Tropical–temperate comparisons in insect seed predation vary between study levels and years

Abstract The biotic interaction hypothesis, which states the species interaction becomes stronger in the tropics, is deeply rooted in classic ecological literature and widely accepted to contribute to the latitudinal gradients of biodiversity. Tests in latitudinal insect–plant interaction have emphasized leaf‐eating insects on a single or a few plant species rather than within an entire community and mixed accumulating evidence, leaving the biotic interaction hypothesis disputed. We aimed to test the hypothesis by quantifying insect seed predation in a pair of tropical and temperate forest communities with similar elevations. We applied a consistent study design to sample predispersal seeds with systematically set seed traps in 2019–2020 and examined internally feeding insects. The intensity of seed predation was measured and further applied to tropical versus temperate comparison at two levels (cross‐species and community‐wide). Our results showed every latitudinal pattern associated with different study levels and years, that is, negative (greater granivory in the tropics in community‐wide comparison in 2020), positive (less granivory in the tropics in community‐wide and cross‐species comparison in 2019), and missing (similar level of granivory in the tropics in cross‐species comparisons in 2020). The cross‐species level analyses ignore differences among species in seed production and weaken or even lose the latitudinal trend detected by community‐wide comparisons. The between‐year discrepancy in tropical–temperate comparisons relates to the highly variable annual seed composition in the temperate forest due to mast seeding of dominant species. Our study highlights that long‐term community‐level researches across biomes are essential to assess the latitudinal biotic interaction hypothesis.

Expanding the research scope to include more species in the community addresses some limitations of single-species studies by increasing taxonomic replication (Anstett et al., 2016). Still, it faces challenges in choosing representative and comparable species sets, thus failing to reach a consensus about the biotic interaction hypothesis (Boyer, 2019;Chen et al., 2017). Furthermore, plant species play different roles in community structure and functioning, and it is questionable to assume that data averaged across selected species can reflect the community-wide spatial pattern (Mottl et al., 2020;Zvereva, Zverev, Usoltsev, & Kozlov, 2020b). The solution is to refocus on the community level where the hypothesis was initially formulated (Coley & Aide, 1991;Coley & Barone, 1996): the percent leaf area of all plant species damaged by insects per year was expected to be greater in tropical forests as a whole, compared to temperate forests. However, the community-wide level tests have been rare, even in insect folivory Rheubottom et al., 2019;Zhang et al., 2011;Zvereva, Zverev, Usoltsev, & Kozlov, 2020b), let alone other types of species interaction.
Compared to leaf herbivory, predispersal seed predation by insects represents a unique but overlooked biotic interaction system (Gripenberg et al., 2019;Janzen, 1971;Lewis & Gripenberg, 2008;Xiao et al., 2017). The predispersal seed predators attack seeds developing on the tree, including members of Lepidoptera, Coleoptera, Diptera, and Hymenoptera (Sallabanks & Courtney, 1992). The insect seed-eaters are engaged in close interactions with host plants and show high specificity (Novotny & Basset, 2005). Insect seed predators have essential ecological functions in inducing host plant fecundity loss and recruitment failure and regulating community composition and local species richness (Gripenberg et al., 2019;Weisser & Siemann, 2008). Furthermore, endophytic insect feeders are directly connected to seed mortality (Lewis & Gripenberg, 2008) and thus reflect the cost of seed predation to the host plants. However, meta-analyses detected absent latitudinal trends in plant and predispersal seed predator interaction (Moles & Westoby, 2003, Zvereva & Kozlov, 2021, warranting further comparison of tropical and temperate seed predation intensity by insects.
We therefore studied seed predation by insects in a pair of tropical and temperate forest communities with similar elevations to test the biotic interaction hypothesis. To our knowledge, this is the first community-level study to compare granivory across biomes.
We specifically asked two questions: (1) Is seed predation intensity higher in tropical forests than in temperate forests? (2) Do the study levels (cross-species vs. community-wide level) influence the results of tropical versus temperate comparisons? 2 | MATERIAL AND ME THODS

| Study sites
We conducted our study in two forest dynamics plots engaged in the Chinese Forest Biodiversity Network (CForBio, Feng et al., 2016) and Forest Global Earth Observatory network (CTFS-ForestGEO, Anderson-Teixeira et al., 2015). Both sites are in protected areas.
The tropical forest plot (XSBN) is located in Xishuangbanna, Yunnan province (21.61°N, 101.57°E), with an area of 20 hm 2 . The elevation ranges from 709 to 869 m, the annual mean temperature is 21.8°C, and the mean annual rainfall is 1493 mm (Lan et al., 2008).
The primary vegetation type is tropical seasonal rainforest, characterized by Parashorea chinensis in the canopy, Sloanea tomentosa, Pometia pinnata, Pittosporopsis kerrii, Garcinia cowa, and Orophea laui in the understory, and Castanopsis echinocarpa near the ridge (Dou et al., 2018). The temperate forest plot (CBS) is located in Mt.
Changbaishan, Jilin Province (42.38°N, 128.08°E), with an area of 25 hm 2 . The elevation ranges from 791 to 809 m, the annual mean temperature is 3.6°C, and the mean annual rainfall is 700 mm (Zhang et al., 2008

| Seed sampling
Seeds were sampled using 150 seed traps with a surface area of ~0.5 m 2 erected in both forest plots following the Barro Colorado Island (BCI) forest dynamic plot (Wright et al., 2005). The seed traps were set systematically and as evenly as possible to ensure community-wide sampling. Each plot was first divided into many 20-× 20-m grids and then separated into four equal subplots by a cross band. The seed traps were finally systematically placed at the grids in the cross band and subplots (Dou et al., 2018;Zhang et al., 2008). The minimum distances between traps were 37 and 20 m for CBS and XSBN, respectively. The distribution maps of seed traps were shown in Dou et al. (2018) for XSBN and Zhang et al. (2008) for CBS.
All seeds and fruits falling into the traps were collected, sorted, and identified to species every 14-15 days. Seeds/fruits collected on the same day were hereafter referred to as a batch of samples.
Seeds/fruits from each trap belonging to one species were then packed, counted, dried, and weighted. We surveyed seeds/fruits across the fruiting seasons of 2019 and 2020, that is, the whole

| Seed predation by insects
We examined all the seeds/fruits through an X-ray machine (Faxitron X-ray Corporation MX-20-DC12, Figure 1) and dissected the suspected seeds under a stereo microscope (OLYMPUS SZ61) to confirm the predation status. Seeds were considered to have been depredated by insects if there were (1) insects (including eggs/larvae/adults) or (2) frass, feeding damage, or entry/exit holes (Chen et al., 2017;Gripenberg et al., 2019;Jeffs et al., 2018). The intensity of seed predation may be derived with different approaches (Basset et al., 2018;Gripenberg et al., 2019;Jeffs et al., 2018), based on either both criteria (1, 2) to reflect the percentage of seed attacked by insects or only the criterion (1) to estimate the potential load of insects on seed mortality. Here, we referred to the former as seed predation rate (proportion of seeds showing signs of seed predator attack) and the latter as incidence of seed predators (proportion of seeds detected with insect predators). We weighted the two metrics by the average seed mass per trap × batch to account for the intraspecific variation in seeds for further analyses. We also retained unweighted metrics to facilitate comparisons with other studies. The incidence of seed predators was calculated for both 2019 and 2020 samples, and the seed predation rate was calculated for the 2020 sample only.
Regarding fruits containing multiple seeds, insect predation was analyzed for each seed separately. The survival of individual seed was our concern in analyses of predispersal seed predation (Chen et al., 2017). On minimal occasions, insect damage was restricted in the flesh and impeded an accurate estimate of seed predation intensity. We counted such a case as one seed predation event, which was unlikely to affect the general results due to the rarity of such cases. Immature seeds were included to minimize the errors in estimating seed predation intensity because premature abscission might be due to insect infection (Jeffs et al., 2018), and insect predation was indeed detected from immature seeds in our samples.

| Simulated historical seed predation rate of Pinus koraiensi in temperate CBS
The dominant species of Korean pine (Pinus koraiensis) in CBS produced seeds in late September and early October, and seed production showed periodic dynamics every 2-3 years (Ji et al., 2002).
After mast seeding in 2018, P. koraiensis produced almost no seeds in 2019 and just a few seeds in 2020. Rodents, birds, and humans F I G U R E 1 The radiograph of Castanopsis megaphylla collected from tropical XSBN in 2020. One intact (bottom left) and two damaged seeds, where a larva (top) and frass (bottom right) were detected, are shown here. The nuts are ~1 cm in length.
further removed some seeds from the trees (Ji et al., 2002;Zhang et al., 2008), and only two seed cones were sampled in 2020. P. koraiensis was underrepresented in our samples, and we used historical data on seed attacks (Heilongjiang Fenglin Nature Reserve, 1975;Mao et al., 2003;Zhang et al., 2020) to simulate the seed predation rate of P. koraiensis by insects before the 1990s when humans began to collect seed cones for sale.
The primary insect predators of P. koraiensis seeds in Mt.
We compared seed predation intensity between tropical and temperate forests at community-wide and cross-species levels.
Community-wide seed predation intensity was defined as the percentage of seed mass of all plant species predated by insects at a given site in a specified period, following Zvereva, Zverev, Usoltsev, and Kozlov (2020b). It was calculated for each batch of samples from both study sites. To test the study sites' effect on ln(1 + √x)transformed seed predation intensity, we implemented a linear model using generalized least squares (gls function in nlme). The errors allowed to be correlated and/or have unequal variances to address temporal autocorrelation and heterogeneity of variance. The batches were coded as the days since 1 January of the year and were included as a covariate. The correlation structure was specified as an autoregressive process of order 1 (correlation = corAR1), indicating the correlation between observations separated by a one-time unit (i.e., 14-15 days) is likely to be more similar than those separated by larger time units (Zuur et al., 2009).
In cross-species levels analyses, we applied the phylogenetic generalized least squares regression (pgls function in caper) to explore the study sites' effect on seed predation intensity while controlling for potential plant phylogenetic nonindependence (Mundry, 2014), which might play a key role in shaping insect herbivory assembly (Turcotte et al., 2014;Weiblen et al., 2006). The phylogenetic signal (measured as λ, lambda = 'ML') was adopted to reflect the extent to which seed predation intensity was statistically related to plant phylogeny (Symonds & Blomberg, 2014). The response in the model was the ln(1 + √x)-transformed seed predation intensity of each species, pooled across all traps × batches per year. The fruit types (fleshy fruit and dry fruit, Yang et al., 2010) and seed mass were possibly confounding variables (Basset et al., 2018;Chen et al., 2017;Gripenberg et al., 2019) and were thus included as covariates. The seeds of gymnosperms, including Pinus koraiensis, Gnetum pendulum, G. montanum, and Abies holophylla, were surrounded by no protective tissue like flesh against insect predators, thus were functionally classified as dry fruits.
Only species with ≥50 seeds were retained in the analyses. The plant phylogeny was constructed with taxonomies against a mega tree as our framework (phylo.maker function in V. PhyloMaker). The plant taxonomy followed Flora of China (Editorial Committee of Flora of China, 1989China, -2013.

| RE SULTS
In the tropical forest plot (XSBN), 15,128 seeds (dry weight: 3958.59 g) from 154 species in 127 genera and 60 families were collected in 2019 (Table 1). The top plant species contributing to seed production was Castanopsis echinocarpa, yielding 14.80% of seed mass (19.43% of seeds). The annual incidence of seed predator was 8.40% (Table 1). A total of 23,589 seeds (dry weight: 4346.75 g) from 159 species in 129 genera and 61 families were collected in 2020.
The top dominant plant species contributing to seed production was Pometia pinnata, yielding 11.67% of seed mass (2.79% of seeds). The annual seed predation rate and incidence of seed predator were 24.73% and 6.72%, respectively (Table 1).
In the temperate forest plot (CBS), 15,614 seeds (dry weight 1093.50 g) from 10 species in four genera and four families were collected in 2019. 79.34% of seeds came from Fraxinus mandshurica, accounting for 52.36% of total seed mass. The annual incidence of seed predator was 17.93% (Table 1). A total of 148,586 seeds (dry weight 3262.18 g) from 14 species in eight genera and eight families were collected in 2020. 91.38% of seeds come from Tilia amurensis, accounting for 74.50% of seed mass. The annual seed predation rate and incidence of seed predator were 5.23% and 2%, respectively (Table 1).
In cross-species comparison, seven species in CBS and 47 species in XSBN with ≥50 seeds were retained in 2019, and nine species in CBS and 56 species in XSBN with ≥50 seeds were retained in 2020. The incidence of seed predator was lower in the tropical XSBN

| DISCUSS ION
Our study has demonstrated that tropical-temperate comparisons in insect seed predation could be higher, lower, or similar at the community-wide (Figures 2 and A1, Tables 1 and 2) or cross-species level (Figures 3 and A2-A4, Tables 1 and 3). We confirm the patterns by two metrics of seed predation intensity in 2020 (i.e., seed predation rate and incidence of seed predators) and by including simulated data of Pinus koraiensis.
Our results show that tropical forests do not consistently suffer greater losses to seed predators than temperate forests, and the study levels and years influence the latitudinal comparisons.
The complex latitudinal patterns thus challenge the generality of the biotic interaction hypothesis, aligning with recent communitywide surveys on granivory, herbivory, and predation on insects (Chen et al., 2017;Mottl et al., 2020;Zvereva, Zverev, Usoltsev, & Kozlov, 2020b Note: Measurements incorporating simulated data of Pinus koraiensis are shown in square brackets. Mean ± SE values are presented. Cross-species measurements unweighted by seed mass, which ignore intraspecific variations in the seed mass, are shown in Table A1.

F I G U R E 2
The tropical-temperate comparisons of the incidence of seed predators at the community-wide level in (

| Study levels
Our findings suggest that the latitudinal pattern of seed predation by insects at the community-wide level does not necessarily match that averaged across individual plant species. The cross-species level analyses ignore differences among species in seed production and weaken (in 2019, between study levels was also demonstrated in insect folivory. Although herbivory values averaged across seven plant species decreased with latitude, community-wide herbivory showed no statistically significant latitudinal pattern (Zvereva, Zverev, Usoltsev, & Kozlov, 2020b). It is thus not reasonable to assume data on the relative losses of seeds to insects collected from single species or averaged across multispecies can be used as a proxy for community-wide seed losses in studies addressing spatial patterns of granivory.
Species-and community-based approaches represent two inter-

| Annual variations
Our findings further suggest that tropical vs. temperate comparisons existed between-year variations at both study levels (Figures 2 and   A1, Figures 3 and A2-A4). We propose the interplay between mast seeding and insect seed predator in temperate CBS as one of the causes. Although mast seeding occurred in both of our study sites (Dou et al., 2018;Qian et al., 2019), it was more pronounced in the temperate CBS, as indicated by a sharp increase of annual seed mass by three times (Table 1). According to the seed predator satiation hypothesis (Janzen, 1971;Kelly, 1994;Silvertown, 1980), plants suffer less seed predation in high-seed years than in low-seed years. The insect predators are easier to satiate where one species dominates local seed production and fewer alternate resources are available, so masting is especially well-developed in temperate forests than in biodiverse tropical forests (Kelly & Sork, 2002;Pearse et al., 2020;Zwolak et al., 2022). Consequently, the seed production and associated seed predation intensity of individual dominant species are more likely to be reflected in the community-wide measurements at higher latitudes. TA B L E 2 Effect of study sites (XSBN vs. CBS) and sampling date (of each batch) on seed predation intensity at the community-wide level

F I G U R E 3
The tropical-temperate comparisons of the incidence of seed predators at the cross-species level in (a) 2019 and (b) 2020. Different patterns are shown here: the incidence of seed predator in tropical XSBN is lower in 2019 but becomes statistically similar in 2020 as against temperate CBS. Violin plots show the distribution of measurements derived from each plant species with ≥50 seeds (CBS: n = 7 and 9 in 2019 and 2020, respectively; XSBN: n = 47 and 56 in 2019 and 2020, respectively), with median and inter quartile range as dots and boxes. Results using seed predation rate as the metric or two metrics unweighted by seed mass are similar and shown in Figures A2 and A3. Results remain similar when including simulated seed predation rate of Pinus koraiensis in CBS, as shown in Figure A4. In this study, the mast seeding of dominant tree species in CBS was so overwhelming that it largely determined community-wide seed predation intensity well below or above tropical measurements. For example, mast seeding of Tilia amurensis in 2020 decreased CBS's community-wide seed predator incidence per batch from 4.10 ± 0.78% (excluding T. amurensis) to 1.84 ± 0.43% (vs. 6.78 ± 0.86% in XSBN, Table 1), whereas mast seeding of Fraxinus mandshurica in 2019 increased CBS's community-wide seed predator incidence per batch from 10.77 ± 3.29% (excluding F. mandshurica) to 17.56 ± 1.07% (vs. 9.73 ± 1.72% in XSBN, Table 1). It is worth noting that the masting seeding, as a population-level adaptive reproductive strategy (Kelly & Sork, 2002), does not necessarily lower seed predation intensity at the community-wide level.
The predator satiation effect of masting is complicated by the intricacies of predator life histories (Zwolak et al., 2022). The life history, mobility, and diet breadth of insect predators and their interactions with other seed predators might shape seed predation patterns and masting dynamics (Bogdziewicz et al., 2021;Kelly, 2021).
However, the insect seed predator guild is poorly studied in our study sites. Further research on insect predators' taxonomies, host specificity, and life histories is needed to comprehend the causes and selective consequences of mast seeding in our study sites.
The annual variations in tropical-temperate comparison (Figures 2 and 3) also indicate that measuring a snapshot of granivory at the community-wide or cross-species level in a specific year might provide opposite and misleading patterns. The among- year variations also exist in leaf herbivory (Adams & Zhang, 2009;Zhang et al., 2011;Zvereva, Zverev, & Kozlov, 2020a). A notable example is insect outbreaks, which are common from a phytocentric perspective and have major top-down effects on plant communities and ecosystems (Carson et al., 2008). We argue the annual variations on some occasions could determine the latitudinal pattern present or absent (Zvereva, Zverev, & Kozlov, 2020a) and change the direction to be positive or negative. Therefore, temporal variation should be considered in testing the biotic interaction hypothesis, and longterm comparisons are recommended (Anstett et al., 2016).

| Interspecific variations
The interspecific variations in seed predation intensity are ubiquitous in the community (Gripenberg, 2018;Jeffs et al., 2018), and understanding the causes is central to studying plant-enemy interaction (Gripenberg et al., 2019). A variety of plant traits and phylogeny have been proposed to influence plant species' susceptibility to predispersal seed predators (Gripenberg et al., 2019). However, neither has received general support at large spatial scales (Basset et al., 2018;Chen & Moles, 2018;Moles & Westoby, 2003), indicating they operate in complex manners or at a fine spatial scale. Our study confirmed such inconsistency and showed phylogenetic signals were moderate to strong (λ = 0.49-0.87, Table 3), seed mass was influential in one study year (i.e., 2020), and fruit types were negligible. The patterns of seed predation intensity across biomes are also likely to be driven primarily by factors not investigated here (such as density-dependent effects mentioned above), and specific research incorporating comprehensive influential factors is further needed.
We statistically accounted for phylogenetic relatedness and two plant traits in the datasets to reduce the variation associated with tropical and temperate plant assemblages. An alternative approach recommended is choosing a subset of plant species from the tropical community that phylogenetically matches the temperate counterpart (Anstett et al., 2016). However, a phylogenetic match does not mean a functional match, as plant functional traits (e.g., seed mass) were not readily incorporated. Seed mass is a central trait in TA B L E 3 Effect of study sites (XSBN vs. CBS), fruit type (fleshy vs. dry fruit), and average seed mass (of each plant species) on seed predation intensity at the cross-species level Note: The phylogenetic generalized linear models were conducted to control plant phylogenetic relatedness. λ was used to measure phylogenetic signal in seed predation intensity, with λ = 0 indicating no phylogenetic signal and λ = 1 indicating a strong phylogenetic signal. Results incorporating simulated data of Pinus koraiensis are shown in square brackets. The p values <.05 are bolded. The results were similar when ignoring intraspecific variations in the seed mass and using unweighted metrics, as shown in Table A2.
seed functional ecology and is associated with seed defense and nutrition (Fenner & Thompson, 2005). Seed mass also exhibits a striking negative latitudinal pattern, declining by 2-3 orders of magnitude between the equator and 60° (Moles et al., 2007;Moles & Westoby, 2003); thus, it is not feasible to choose a subset of plant species with comparable seed mass with temperate counterpart.
Moreover, phylogenetically selected species usually account for a tiny proportion of plant assemblage in the biodiverse tropical forest, making representativeness insufficient.

| CON CLUS ION
We applied a consistent community-wide sampling design in a pair of tropical and temperate forests in mainland East Asia to study predispersal seed predation by insects. Our results reveal every latitudinal pattern associated with different study levels and years, to some extent reflecting the controversy surrounding the biotic interaction hypothesis in a single study. The latitudinal pattern of species interaction could be far more complicated due to interspecific, temporal, and other variations (Anstett et al., 2014;Coley & Barone, 1996). Like in testing the density-dependent effects (Bogdziewicz et al., 2021;Cannon et al., 2021), incorporating and addressing high and pervasive variation in natural systems is vital in developing a realistic assessment of the latitudinal biotic interaction hypothesis.

ACK N OWLED G M ENTS
We thank the staff of the National Forest Ecosystem Research Station

CO N FLI C T O F I NTE R E S T
None declared. Note: The phylogenetic generalized linear models were conducted to control plant phylogenetic relatedness. λ was used to measure phylogenetic signal in seed predation intensity, with λ = 0 indicating no phylogenetic signal and λ = 1 indicating a strong phylogenetic signal. Results incorporating simulated data of Pinus koraiensis are shown in square brackets. The p values < .05 are bolded. Results with weighted seed predation intensity, which incorporated the intraspecific variations in seed mass, are shown in Table 3.

F I G U R E A 2
The tropical-temperate comparisons of seed predation rate at the cross-species level in 2020. Tropical XSBN is similar in seed predator rate to temperate CBS, regardless of the metric being weighted by seed mass (a) or not (b). Violin plots show the distribution of measurements derived from each plant species with ≥50 seeds (CBS: n = 9; XSBN: n = 56), with median and inter quartile range as dots and boxes. Results using the incidence of seed predators as the metric are similar and shown in

F I G U R E A 4
The tropical-temperate comparisons of seed predation rate at the cross-species level in 2020, including simulated data of Pinus koraiensis in CBS. Tropical XSBN is similar in seed predator rate to temperate CBS, regardless of the metric being weighted by seed mass (a) or not (b). Violin plots show the distribution of measurements derived from each plant species with ≥50 seeds (CBS: n = 10; XSBN: n = 56), with median and inter quartile range as dots and boxes. Results using real data are similar and shown in Figures 3, A2, and A3.

F I G U R E A 1
The tropical-temperate comparisons of seed predation rate at the community-wide level in 2020. The seed predation rate in tropical XSBN is higher than that in temperate CBS (a). The result remains similar when replacing (thus increasing) the seed predation rates of late September and early October batches in CBS with the simulated data of Pinus koraiensis (b). Violin plots show the distribution of measurements derived from each batch of seed samples (CBS: n = 13; XSBN: n = 27), with median and inter quartile range as dots and boxes. Results using the incidence of seed predators as the metric are similar and shown in Figure 2.